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ABSTRACT 

The accreting millisecond X-ray pulsar SAX J1808.4-3658 shows peculiar low luminosity states known 
as “reflares” after the end of the main outburst. During this phase the X-ray luminosity of the source varies 
by up to three orders of magnitude in less than 1-2 days. The lowest X-ray luminosity observed reaches a 
value of ^lO^^ergs^*, only a factor of a few brighter than its typical quiescent level. We investigate the 2008 
and 2005 reflaring state of SAX J1808.4-3658 to determine whether there is any evidence for a change in the 
accretion flow with respect to the main outburst. We perform a multiwavelength photometric and spectral study 
of the 2005 and 2008 reflares with data collected during an observational campaign covering the near-infrared, 
optical, ultra-violet and X-ray band. We And that the NIR/optical/UV emission, expected to come from the 
outer accretion disk shows variations in luminosity which are 1-2 orders of magnitude shallower than in X- 
rays. The X-ray spectral state observed during the reflares does not change substantially with X-ray luminosity 
indicating a rather stable configuration of the accretion flow. We investigate the most likely configuration of 
the innermost regions of the accretion flow and we infer an accretion disk truncated at or near the co-rotation 
radius. We interpret these findings as due to either a strong outflow (due to a propeller effect) or a trapped disk 
(with limited/no outflow) in the inner regions of the accretion flow. 

Subject headings: neutron stars - pulsars - X-ray binary - accretion, accretion disks - X-rays/pulsars: individ¬ 
ual (SAX J1808.43658) 


1. INTRODUCTION 

SAX J180 8.4-3658 is a binary X-ray transient discov¬ 
ered in 1996 (lin_TZand_eLalJ[I998|) a nd located at a distance 
of 2.5 -3.5 kpc (iin ’t Zand et al.lll998tlGallowav & Cummingl 
|2006|) . The binary is c omposed by a 401 Hz accreting mil- 
lisecond X-ray pulsar (IWiinands & van der KlisI Il998t see 
iPatruno & Watts! 120121 for a r eview) and a semi-degenerate 
companion o f mass ^0.07Mp) (Bildsten & Chakrabartv 20M1 : 
iDelove et akl 120081). The system has a n orbital pe¬ 
riod of 2 hr (ICh^abartv & MorganI Il998h and the ac¬ 
creting neutron star has an inferred magnetic held of 
^10^ G (iPsaltis & Chakrabartvl [19991 iDi Salvo & Burderil 

120031 [Hartman et al.ll2008l: IPatruno et alJl201^ 

The seven outbursts of SAX J1808.4—3658 observed so 
far (in 1996, 1998, 2000, 2002, 2005, 2008 and 2011) have 
a recurrence time of ^^2-3.5 years and shows a main out¬ 
burst with a fast rise followe d by a slow decay an d a rapid 
drop (see e.g.. Figure 2 in iHartman et akl l2008h . After 
the end of the main outburst, lasting for approximately one 
month, there is a long “reflaring” taiQ at low X-ray lumi- 
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nosities (10 ergs ) lasting for several 

e.g.JWiinands et al.(l2001t iRevnivtsevI 12003 

tens of days (^see 

lWiinandsll2003l; 

IHartman et al.1120081 IPatruno et al.112009bl 

2012^. Such re- 


flares manifest themselves as smooth variations in luminos¬ 
ity that appear cyclically every few days. Furthermore, on 
one occasion (during the 20 05 outburst) a very lo w luminos¬ 
ity phase has been reported (ICampana et aklUOP^ . This later 
phase was called “metastable” since it occurred at an X-ray lu¬ 
minosity of 1.5 — 5 X 10^^ergs^' (0.5-10 keV), a value com¬ 
patible with the lowest X-ray luminosities reached d uring the 
reflares (IPatruno et al.ll200^ IWijnands et al.ll2()()TI) and only 
a factor of a few brighter than the typical quiescent luminosi¬ 
ties. It is currently unclear whether this phase is really differ¬ 
ent from the reflaring phase or whether the sparse Swift mon¬ 
itoring observed only the low luminosity phase of the reflares 
themselves. After this possible metastable state the source 
turned back to quiescence, when the 0.5-10 keV X-ray lumi- 
nosity reached a value of 5 x 10^'ergs^' (ICampana & Stell^ 
12001 iHeinke’eT^IIOO^ . 

The typical variation on luminosity observed during the re¬ 
flares of SAX J1808.4—3658 spans approximately one or- 
derof magnitude (10^^— lO^^ergs^', see e.g.. IPatruno et alJ 
|20093). However, in at least three occasions (during the 2000, 
2005 and 2008 outbursts) the reflares showed a dramatic de¬ 
crease in luminosity with an excursion of thr ee orders of mag¬ 
nitude on a tirnescale shorter than ^2 day s (IWijnandsll2003l : 
ICampana et alJ[2008t IPatruno et al.ll2009bl) . The faintest lu¬ 
minosities reach a value of '^3 x lO^^ergs^* (for a distanc43 
of 3.5 kpc) which is slightly higher than the typical quies- 

[200^. _ 

^ IWiinandl <200.311 reports a luminosity of 1.7 x 10^’ergs^* foradistance 
of 2.5 kpc. Rescaling to 3.5 kpc - the distance adopted in this work - gives 
3.3 X 10^^ergs“^ 
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cent one (5-8xl0^*erg s~^ at 3.5 kpc. iHeinke et al.l 120091: 
ICampana & Stella! I2004ft . Such low luminosities were ob¬ 
served thanks to the sensitivity of the XMM-Newton and Swift 
observatories and we cannot exclude that something similar 
might have happened in all the other outbursts recorded so 
far (which were observed, in X-rays, exclusively by the Rossi 
X-ray Timing Explorer, RXTE). 

Reflares (also known in the literature as “rebrighten- 
ings”, “echo-outbursts”, “mini-outbursts” and “flaring-tail”) 
have been observed in a number of different systems, 
including dwarf nova e (WZ Sae, iPatterson et alj |2002|: 
EG CnC. lOsaki et alJ 1200 U ER U MA Robertson et alJ 
Il995t IPatterson et alJ 1201 3L AL Com [Howell et al. 19961 
UZ Boom. iKuulkers et al.l 1199'^ BK Lvn jPatterson et M 
1201 3h . neutron star s (e.g., KS 1731 -260, iSimonI 120101 
SAX J1750.8- 2900 lAllen etal.l ImsI and possibly IGR 
J00291H-5934, iLewis et al.l |20 1C ) and black hole low 
mass X-ray binaries (e. g., XTE J1650-500. iTomsick et ^ 
[2001 GRO .T0422-H32. kuulkers et al.1119961: GRS 1009- 
45. iBailvn & Oroszl 119951: XTE J1859-H226. iZurite et alJ 
I2002lk The ubiquity of the phenomenon among different 
types of accreting compact objects suggests that reflares are 
related to specific properties of the accretion disk itself rather 
than properties of the compact object. 

Reflares are not a feature that can be easily justified within 
th e disk instability model (DIM; see the extended reviews 
oflSm^ll984t iMenou et al1l2000t [Pubus et al]l2001l : iLasotal 
120011: iKotko et al.l 120121) . Reflares are observed in numeri¬ 
cal DIM simulations and they are created by back and forth 
propagation of cold and hot waves that re-combine and ion¬ 
ize the accretion disk cyclically. However, these reflares are 
considered a “deficiency” of the DIM model rather than a fea¬ 
ture, since they cannot reproduce almost any of the proper¬ 
ties observed in dwarf no vae and soft X-ray transients (iLasotal 
l2001l:[Dubus et al.l2001ll . Eurthermore reflares require a large 
reservoir of matter to be still in the disk once the main out¬ 
burst is over, whereas the DIM predicts that the accretion disk 
shoul d be almost devoid of matter toward s the end of an out¬ 
burst (lArai et al.ll2009l : [Dubus et al1l2001[) . 

An interesting aspect of the reflares seen in 
SAX J1808.4—3658 is that the physical conditions in 
the inner regions of the accretion disk are strongly affected by 
the presence of its magnetosphere, since accretion-powered 
pulsations are observed throughout the outbursts, including 
large portions of the reflares in 2000, 2002, 2005 and 2008 
(at least at luminosities ~ 10^^ergs~ ^, where RXTE has 
suffic ient sensitivity to detect them, see Hartman et al.ll2008l 
[ 200 ^ . 


During the 2008 outburst, SAX J1808.4—3658 had been 
monitored with a multi-wavelength campaign that involved 
the RXTE, the Swift X-Ray Telescope (XRT) and Optical and 
Ultra-Violet Telescope (UVOT) and by ground-based opti- 
cal/NIR observations in the I and H bands taken with the 
1.3-m SMARTS telescope. The fact that SAX J1808.4-3658 
is an accreting millisecond pulsar with a dynamically impor¬ 
tant magnetosphere plus the unique multiwavelength cover¬ 
age is key to simultaneously probe different regions of the 
accretion disk/flow. In this paper we report the results of the 
multi-wavelength campaign on the 2008 outburst with a par¬ 
ticular focus on its reflares. We also analyze the 5w//f/UVOT 
and XRT data collected during the 2005 reflaring phase and 
“metastable” state and compare them with th e 2008 outburst 
and with the 2005 analysis already reported in ICampana et alJ 


( 120081) . We investigate whether during the reflaring phase 
our observations show any evidence for anomalies in the be¬ 
haviour of the accretion flow and whether there is any evi¬ 
dence in support of any of the current models that try to ex¬ 
plain the reflares phenomenon. 

2. OBSERVATIONS 

We used data collected with four different detec¬ 
tors/telescopes: the optical/NIR ANDICAM on the 1.3- 
m S MARTS telescope (oper ated by the SMARTS Consor- 
tium;ISubasavage et al.l20i^ . the 5 w/ff/UVOT (iRoming et alJ 
120051) , the SwiftlXKY (0.3-10 keV: iBurrows et alJl2005[) ^ 
the Proportional C ounter Array (PCA) aboard RXTE (2-60 
keV: iJahoda et alJ 120061) . In the following we describe the 
data reduction procedure for each waveband. 

2.1. X-Ray Data 

We analyzed 20 targeted Swift observations taken between 
September 18 and November 8, 2008 and 21 taken between 
June 17 and 28 October 2005 (see Table [TJ. Ten observa¬ 
tions of 2008 and six of 2005 were taken during the main out¬ 
burst, whereas the remaining monitored the reflaring phase. 
The XRT was operated in either photon-counting (PC) mode 
(2.5073-s time resolution) or in window-timing (WT) mode 
(1.7675-ms time resolution). 

We extracted Swift count rates, spectra, and spectral re¬ 
sponse files for each individual ObsID using the online 
5w/ff/ XRT data products generatoiF^ (lEvans et al.l200’^l2009L 
120141) . The spectra were extracted (and fitted) in the 0.3- 
10 keV band using the default event grades (0-12 for PC 
mode and 0-2 for WT mode). The spectra were grouped with 
the GRPPHA task in HEASOET (ver. 6.16) and fitted with 
XSPEC (ver. 12.8.2). Spectra containing >150 counts were 
in general grouped to a minimum of 20 counts per bin and 
fitted with the statistic, whereas spectra with fewer counts 
were grouped to a minimum of 1 count per bin and fitted with 
the W statistic (a modified version of the C statistic that al¬ 
lows background to be taken into account). Eor a few spectra 
with the largest number of counts we used a higher group¬ 
ing minimum than 20 to avoid oversampling the spectral res¬ 
olution of the detecter by a large factor. Spectra with <15 
counts were not fitted; instead the count rates were converted 
to fluxes assuming an absorbed power-law model. Eor Ob- 
sIDs where both PC and WT mode data exist we fitted spectra 
from both modes simultaneously (with all parameters tied) if 
the numbers of counts in the spectra were of the same order 
of magnitude, but otherwise only used data from one or the 
other mode (which was usually the case). We also calculated 
a hardness ratio for each ObsID based on the PC mode data 
only (where available); this was defined as the net counts in 
the 2-10 keV band divided by those in the 0.3-2 keV band. 

Most of the data recorded in 2005 and 2 008 by RXTE 
was al read y presented iniPatruno et aT] (I2009bl) . iPatruno et akl 
(120121) and iHartman et al.l (120091) andwe refer to those works 
for in-depth details. Here we report a summary of that data 
analysis which is relevant for the results presented in this pa¬ 
per. Eor the RXTE data, we extracted the 2-16 keV energy 
band flux from the Standard2 mode data (16-s time resolution) 
collected by the PCA instrument during the 2005 and 2008 
outbursts. The background counts were calculated with the 
ETOOL pcabackest and were subtracted from the lightcurve 

http://www.swift.ac.uk/user_objects/ 
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along with dead-time corrections. The energy-channel con¬ 
version was done by using the pca_e2c_e05v04 table pro¬ 
vided by the RXTE team. 

2.2. Space-borne UV/Optical Observations 

The 5w//t/UVOT operated in imaging mode during all ob¬ 
servations. Most of the observations (during the 2008 out¬ 
burst) were taken with six filters (v, b, u, wl, m2 and w2) 
wheres the last 4 exposures used the u filter only. The 2005 re¬ 
flare were always observed with the v filter. We extracted the 
source photons from a circular region with a radius of 2.5 arc- 
seconds. The background was extracted from a circular region 
with radius of 10 arcseconds far from bright sources. Before 
proceeding further, we manually inspected all UVOT images 
and found that some snapshots were affected by poor track¬ 
ing. In particular some snapshots within the sequences corre¬ 
sponding to ObsIDs 30034036, 37 & 38 were affected and we 
only used those snapshots that were unaffected. Since the data 
was taken in image mode, bad snapshots themselves could not 
be corrected. We used the tool uvotimsum to add all UVOT 
exposures (snapshots) present in each FITS file into a single 
high signal-to-noise image and uvotsource to d etermine the 
optica l/UV magnitudes using the Vega system (iPoole et alJ 
l2008l) . We then proceeded to de-redden the magnitudes by 
choosing a value for the neutral hydrogen absorption column 
of Nh = (1.4±0 2) X lO^^cm"^ dP^no et alJl2009i) . We 
used the relation (iGiiver & Ozelll200^ : 

V/r = (6.86 ±0.27) X 102'£(B-V) (1) 

to obtain E{B — V) = 0.204 ± 0.030. The re ddening fo r the 
six different filters is calculated according to iPeil (Il992h (see 
Tablelali. 

2.3. Ground-Based Optical and Near-Infrared Observations 

The ground-based NIR (/ and //-band) observations 
were made using the 1.3m SMARTS telescope at the 
Cerro Tololo Inter-American Observatory (CTIO) in Chile. 
SAX J1808.4—3658 was observed with a roughly daily 
cadence (weather permitting) between September 22 and 
November 11, 2008. The /-band data were anal yzed with the 
standard IRAF optical data reduction pipelines (IBuxton et alJ 
1201 2|) whereas for the NIR data multiple dithered frames 
were taken and then flat-fielded, sky subtracted, aligned, and 
average-combined using an in-house IRAF script. Three stars 
nearby SAX J1808.4-3658 and within the SMARTS FOV 
were used as references and their average magnitudes were 
used as a basis for differential photometry with respect to 
SAX J1808.4—3658. The H-band magnitudes of these refer¬ 
ence stars were taken from the 2MASS point source catalog. 
For I -band we used t he ref erence stars and their magnitudes 
from iGreenhill et ^ (l2006h . and compared these with the in¬ 
strumental magnitudes obtained from SMARTS frames. We 
used the following zero-point fluxe s to convert from m agni- 
tude to flux-d ensity: /n = 2416 Jv (iBessell et al.lll99^ . and 
Hq = 980 Jy (IFrogel et alJlf^ lElias et alJll98a . We also 
de-reddened the magnitudes following the same procedure 
outlined above (see Table|2l. 

2.4. Contamination from Field Stars 

Since the point-spread-function of the UVOT telescope 
is rather poor (2.5 arcsec at 350 nm) and the held of 
SAX J1808.4—3658 is crowded we need to evaluate the 
amount of contaminating light from nearby stars that fall 



Figure 1. UVOT image of the FoV around SAX J1808.4-3658 taken with 
the V filter taken around MJD 53627.17 (ObsID 30034018, Sept. 14, 2005). 
The dashed red circle highlights the 2.5” region around a bright star whose 
light “leaks” into the extraction region used for SAX J1808.4—3658 (white 
circle). 


within our circular aperture extraction region. By comparing 
the UVOT fleld-of-view (FoV ) of SAX J1808.4-3658 with 
Figure 2 in IWang et akl (l2009h we found that part of the light 
emitted by a rather bright star leaks into the extraction region 
ofSAX J1808.4-3658 (see Figure^. 

To verify how much the contamination is affecting our mea¬ 
surements we first selected three observations: one taken in 
2005 (ObsID 30034018), one in 2008 (ObsID 30034040) and 
one during quiescence in August 2014 (ObsID 30034062). 
We chose the 2005 observation because the high signal-to- 
noise of the image allows a precise determination of the cen¬ 
troid of the contaminating object. The observation taken in 
2008 is used to cover all six UVOT Alters whereas the 2014 
one is used to measure the UV flux during quiescence. We 
find that the contaminating object has the following magni¬ 
tudes: ~19.1 in V, > 19.45 in b, ^19.8 in u and > 21 in wl. 

Since the minimum magnitude detected (for the extraction 
region centered around SAX J1808.4—3658) is: ^ 20 in v, 
19.4 in b, '^18.7 in u and 19.3 in wl we can confidently ex¬ 
clude that the contaminating source has an important effect in 
any but the v band. 


3. RESULTS 

3.1. The Multi-Wavelength Lightcurve 

The multi-wavelength lightcurve of the 2008 outburst is 
shown in Figure |2l The reflares are observed at all wave¬ 
lengths (with the exception of the first reflare, not monitored 
by Swift) with a typical cadence of about one observation 
every 2 days for Swift/XRI and UVOT, every 1-2 days in 
NIR with SMARTS and multiple observations per day with 
RXTEfPCA. The main outburst is fully covered with a high 
cadence monitoring with RXTE and observed by SMARTS for 
only two days close to the peak of the outburst and by Swift 
for 10 consecutive days during the slow decay. 

The first 5w(/f/XRT reflare observation is a non-detection 
and corresponds to an upper limit (95% confidence level) of 
less than 5-10 times the value detected in quiescence {Lq ~ 5- 
8 X 1 0^ * erg s^ \ ICampana & Steiyi2004t iHeinke et al.ll200'^ 
120091) . The timescale of the reflares is about 7-12 days, 
which is substantially longer than the fast decay timescale 
(^3 days; for a detailed discussion of the slo w/fast decay in 
SAXJ1808.4-3658 see lHartman et alJ2008l) . By comparing 
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Table 1 

Observations of SAX J1808.4—3658 from 2005 to 2008 


Year 

Instmment 

Program IDs 

Obs-IDs 

Data Range [MJD] 

2005 

RXTEf?CA 

91418,91056 

91056-01-*, 91418-=^-*-* 

53522.8 to 53581.4 

2008 

RXTEf?CA 

93027 

93027-*-*-* 

54731.9 to 54775.2 

2005 

SwifilXKV (PC) - Swiftf\}\OT (y) 

30034, 30075 

30034001 to 30034025 
30075001,30075020 

53538.0 to 53671.0 

2008 

5vy(/r/XRT (PCAVT) - 5vy////UVOT {v,b,u,w},m2,w2) 

30034 

30034026 to 30034044 

54734.0 to 54778.0 

2008 

^MARraCTIO-l.Sm (H, I) 

N/A 

N/A 

54732.1 to 54781.0 


Table 2 

Reddening coefficient for each of the eight optical/UV filters on UVOT and 
SMARTS (E{B -V)= 0.204 ± 0.030). 


Filter 

Reddening {A^ ± a^) 

H 

0.11 ±0.02 

I 

0.37 ±0.06 

V 

0.62 ±0.09 

b 

0.79±0.12 

u 

1.00 ±0.15 

wl 

1.38 ±0.20 

m2 

1.96 ±0.29 

w2 

1.64 ±0.24 


the behaviour of the I and H magnitudes with that of the X- 
rays it appears that the NIR luminosity is higher than what 
would be expected when extrapolating the behavior seen dur¬ 
ing the slow decay. Indeed the two NIR observations occur¬ 
ring right after the beginning of the fast decay have an H and 
I magnitude higher than what would correspond to the X-ray 
flux detected on the closest dates (see Figure |2]i. Therefore 
it is possible that reflares start earlier in NIR by at least 1.5 
days. The beginning of the other two reflares is ill constrained 
since the sampling is not sufficiently dense in X-rays, but the 
observations are still compatible with an earlier response of 
the NIR. 

The last detection of the reflaring phase occurs at 
MJD 54781 (/ band) at a level compatible with what observed 
during the reflares minima. The source is not detected in the H 
band and no observations are performed during that time with 
the SwiftfXKI and UVOT. It is therefore difficult to determine 
whether we have detected the beginning of a new reflare, the 
beginning of a “metas table” state similar to what reported by 
iCampana et al.l (120081) at the end of the 2005 outburst or some¬ 
thing different. The I band magnitude detected at MJD 54781 
is ^2.5 mag. brighter than the qui escent I detected in the 
1998 outburst (ICampana et al.l 1^041) so that we can at least 
confidently conclude that some activity is still ongoing at the 
time of the observation. We also inspected the publicly avail¬ 
able lightcurves taken with the All Sky Monitor aboard RXTE 
and the Burst Alert Telescope on Swift. We found a few data 
points close to the 3(7 detection level, but none of them is sta¬ 
tistically significant given the number of data points (trials) 
considered (i.e., 3248 data points). The constraints on the lu¬ 
minosity placed by these two instruments are not particularly 
strong since they correspond to luminosities of the order of a 
few times lO^^ergs^', which are values reached only during 
the brightest portions of the reflares. 

The V band magnitudes (not dereddened) of the 2005 out¬ 
burst are shown in Figure[3along with the 2-16 keV and 0.3- 
10 keV X-ray flux. Several observations (highlighted with 
blue circles in the figure) show a lack of correlation between 
the V magnitudes and the X-ray flux. 


3.2. X-Ray-Optical/UV/NIR correlation 

By looking at Figure|2]we can see how the X-ray luminos¬ 
ity variations co-vary with the the NIR, optical and UV lumi¬ 
nosity. However, the excursion seen in X-rays spans at least 
2 orders of magnitude, whereas the NIR/optical/UV changes 
by approximately 3 magnitudes, which translates into a flux 
variation of a factor r:! 10-20. Therefore, we do not expect 
a steep correlation between X-ray flux and NIR/optical/UV 
magnitudes. Furthermore, on several occasions the X-rays 
and NIR/opticaFUV do not correlate and seem to be varying 
independently (see e.g., circled points in Figure [3] and data 
poi nts aro und MJD 54760 in Figure |2]and discussion in Sec¬ 
tion |4TT]l. 

In this section we estimate the correlation between these 
quantities. Since the variations in luminosity can happen on 
very short timescales, we decided to consider only the optical, 
UV and X-ray emission recorded with the UVOT and the XRT 
and to exclude the ground-based optical/NIR data since they 
are not taken simultaneously with the X-ray data, but have 
a typical offsets of ^0.5-1 day. The advantage of using the 
UVOT and XRT data is that these telescopes operate strictly 
simultaneously (although the exposure times are not necessar¬ 
ily the same). One exception is represented by one ground- 
based I and H band observation from MJD 54758 taken 
within 1 hour of a Swift pointing. In this case we decided 
to consider the ground-based data as (quasi)simultaneous. 

In Figure |4] we plot the de-reddened optical to X-ray (un¬ 
absorbed) luminosity relation for the Swift observations of 
the 2005-2008 outbursts plus additional data referring to 
other neutron star low mass X-ray binaries where these two 
quantities are measured (and found to be correlated, see 
iRussell et al.ll2006l 120071 for details). We used the 2-10 keV 
energy ba nd rather than the 0.3 -10 keV in order to match the 
analysis of iRussell et al.l (l2006h done for the other neutron star 
LMXBs used in Figure |4] Some of the 2008 data lie slightly 
above the correlation with a slope that is somewhat flatter than 
recorded for other neutron star LMXBs. The offset between 
our data points and the other archival data is not surprising 
since the latter points are taken with the BVRI Alters, whereas 
our data points are mainly composed by u and UV data and 
only very few data points are taken with the /, b and v Alters. 
Indeed such effect is not observed in 2005 when the UVOT 
used only the v Alter. Similar (instrument al) offsets have been 
obser ved already in other sources (see e.g. Ivan der Horst et af] 

ISoil. 

The flatter slope instead needs an explanation and to quan¬ 
tify this, we first fitted a simple power law relation FjjyjQpj = 

kF^, where Fuyiopt is the flux in one of the optical/UV 
bands, Fx is the 2-10 keV X- ray flux, /3 is the slope an d k 
is a constant of proportionality. [Russell et all (120061120071) es¬ 
timated the coefficient j3 for neutron stars (assuming a radia- 
tively efficient accretion flow. Lx M) and black holes (with 
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Figure 2. Multi-wavelength lightcurve of the 2008 outburst. The X-ray flux and the NIR/optical/UV magnitudes are not corrected for absorption/reddening but 
refer to the observed values. Top panel: the observed X-ray lightcurve recorded with RXTE (2-16 keV; filled black circles) and ^wz/r/XRT (0.3-10 keV; open 
gray circles). The vertical dotted line marks the approximate beginning of the rapid drop. The dashed vertical lines identify the most likely limits for the start and 
end of each reflare. The down-arrow marks a non-detection. Middle panel: optical/NIR observations with the H (central A = 16500A ) and I (central A = 7980A) 
filters. Bottom panel: optical and UV observations carried with the UVOT (for the central A of each filter we refer to Table [sj. The non-detections are not 
reported. 


Table 3 

Correlation coefficient between the UV/optical and X-ray fluxes 


UVOT/Band 

^central (-^) 

P 

w2 

1928 

0.28 ±0.05 

m2 

2246 

0.24 ±0.06 

wl 

2600 

0.22 ±0.05 

u 

3465 

0.23 ±0.05 

b 

4392 

0.16 ±0.04 

'■’2008 

5468 

0.03 ±0.05 

''2005 

5468 

0.05 ±0.21 

/ 

7980 

0.11±0.13 


a radiatively inefficient accretion flow, Lx hfi). Since we 
are using optical/UV wavebands, the expected range of p for 
a viscously heated disk is 0.5 <P< 0.67 for neutron s tars 
and 0.25 < P < 0.33 for black holes (iRussell et al.ll2006h . In 
Table|3]and Figure|5]we report the results of our fit. 


All bands, with the exception of the v and I bands, show 
a flux that is clearly correlated with the X-ray flux. How¬ 
ever, the value of the correlation coefficient j3 falls within 
the 0.15-0.25 range (again, excluding the v and I bands). 
These values are significantly s maller than those seen in 
other neutron sty LM XBs (e.g.. I Armas Padilla et alJ 1201^ 
iMaitra & Bailvnll2008h . The values of p also increase with 
decreasi ng wavelength, as found in viscously heated disks 
(see e.g. IRussell et al.ll200^ lArmas Padilla et al.ll2013h . For 
irradiated disks or for jet-dominated optical emission, the ex¬ 
pected values of P are also higher than observed. 

3.3. Rayleigh-Jeans limit 

The shallow slope of the Fx vs. Fgpfjuy relation could also 
be explained if the optical/UV emission falls in the Rayleigh- 
Jeans (RJ) limit of the multicolor-blackbody accretion disk 
spectrum. To test this hypothesis we selected five observa¬ 
tions with the source detected in at least four UVOT filters. 
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Figure 3. Top Panel: v band lightcurve (SwiftfUVOT) of SAX J1808.4-3658 
during the 2005 outburst reflare. The magnitudes are not dereddened. Bottom 
Panel: 2005 outburst X-ray lightcurve observed with RXTEfPCA (2-16 keV) 
and SwiftlXm (0.3-10 keV). 

Of these observations, four were taken during the reflares 
whereas a “control” group was chosen from the main outburst. 
We then fitted ^nd checked whether atal as ex¬ 

pected in the RJ limit. In one observation we included also the 
ground-based optical/NIR observations since these data were 
recorded within 1 hour from the UVOT and XRT data. We 
also fitted a correlation with a forced to be emal to 2 as ex¬ 
pected in the RJ limit. In Table |4] and Figure we report the 
observation selected and the results of our fit. 

None of the observations show an a compatible with being 
2. Three out of five observations cannot be fitted with a simple 
linear relation (i.e., unaccmtable However, by looking at 
the panel “a)” of Figure 0 which refers to the observation 
taken during the main outburst, it appears that we might still 
be in the RJ limit for all optical/UV filters if the location of the 
data points is somehow affected by unaccounted systematics 
(see discussion below). Panels “d)” and “e)” are compatible 
with the RJ limit if the peak of the black-body emission falls 
around the m2 and w2 band. In the panel “b)” and “c)” there 
seems to be stronger deviations from the RJ limit, with a be¬ 
ing closer to values of 0-1. These two observations might 
be explained if the blackbody peak towards lower energies. 
However, the optical/UV luminosity is slightly lower than that 
observed in panels “d)” and “e)”, which is the opposite effect 
that one would expect. 

An obvious source of systematic uncertainty is the value 
of E{B — V) (dependent on the value of Nr) which we have 
assumed to be constant throughout the outburst. To verify 
whether this is indeed affecting our results we fitted our obser¬ 
vations with the model reddening*diskbb via the software 
XSPEC (v. 12.8.1) by leaving E{B — V) asa free parameter and 
fixing Aff = 1-4 X 10 ^* cm^^, 7i„ =0.1 keV and i = 50° (e.g., 
iPatruno et al.ll200^ . In this way we also take properly into 
account the value of reddening as a function of wavelength 
over the entire filter bandwidth. The first two observations 
have best-fits that are statistically unacceptable with null hy¬ 
pothesis probabilities <C 0.01. Changing the temperature of 



31 32 33 34 35 36 37 38 

log(X-Ray Luminosity [2-10 keV]) 

Figure 4. Top Panel: NIR-optical-UV/X-ray correlation diagram. The 
green p oints show the correlati on for several different neutron star 
LMXBs IRussell et al.ll2006L 1200711 . A few points (circled in red) refer to 
(archival) SAX J1808.4-3658 data taken during quiescence or during its 1998 
outburst. The 2008 Swift observations reported in this paper are shown with 
different symbols/color depending on the specific filter used. Bottom Panel: 
Same as top panel but for the 2005 outburst. 

the diskbb model does not change the results in any significant 
way. The third and fourth observations are compatible with 
a value of E{B — V) ~ 0.01-0.16, which app ears to be only 
marg inally compatible with the galactic Nr (iKalberla et al.l 
120051) . The last two observations give a range of reddening 
which is unconstrained. Therefore we conclude that the data 
cannot be explained with a simple multicolor-disk blackbody 
and that there is no evidence for the data to be close to the 
RJ limit and therefore the low values of j3 need a different 
explanation. 

Finally, we show that indeed the NIR/optical/UV emission 
is always larger (by a factor ^10) than would be predicted 
assuming that the X-ray emission tracks the accretion rate so 
that Lx — GMM/Rt, and that the outer disk is a standard 
Shakura-Sunyaev a-disk. This is seen in Figure |7] which 
plots the observed NIR/optical/UV flux of the 2008 outburst 
between MJD 54757 and 54770 and the predicted optical 
disk emission based on the X-ray luminosity (although note 
that this is based on the 2-10 keV X-ray emission, which 
could underpredict the bolo metric luminosity by a factor of 
up to 3 lin’t Zand et al.ll2007h . The optical emission is brighter 
than would be predicted by X-rays, by a factor of a few (but 
it would track more closely if the X-ray bolometric flux is 
higher). 
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Table 4 

Observations Used to test for the Rayleigh-Jeans limit 


Panel 

Phase 

Obsid 

MJD 

Filters 

a 

Z^/dof 

a 

Main Outburst 

30034031 

54738.4 

v,b,u,wl,m2,w2 

1.4±0.3 

37.7/4 

b 

Reflai'e 

30034036 

54758.0 

H, I, b,u,wl,m2,w2 

0.93 ±0.11 

30.5/5 

c 

Reflare 

30034037 

54760.3 

v,b,u,wl,m2,w2 

0.48 ±0.18 

7.37/4 

d 

Reflare 

30034038 

54762.1 

b,u,wl,w2 

0.58 ±0.28 

2.9/2 

e 

Reflare 

30034039 

54764.8 

b, u,wl,m2,w2 

0.32 ±0.50 

10.80/3 


0.7 


0.6 

0.5 

0.4 


=3- 0.3 


0.2 

0.1 

0 

- 0.1 


Figure 5. Correlation coefficient j3 calculated for seven different optical/UV 
bands for the 2008 outburst. The P value for the v band of the 2005 outburst 
is unconstrained and is not displayed. The two colored bands indicate the 
theoretical P values expected for black hole LMXBs (in green) and neutron 
stars LMXBs (in yellow). The difference in the location of the two bands 
reflects the different accretion efficiency (« M for neutron stars and « for 
black holes). Most data points of SAX J1808.4—3658 fall within the black 
holes band. The value of p is also consistent with increasing towards shorter 
wavelengths, a typical behaviour expected in viscous heated accretion discs 
rather than irradiated ones. 

3.4. X-Ray Spectra and Colors 

We have analyzed the X-ray spectra and performed a color 
analysis of SAX J1808.4-3658 during the 2008 and 2005 re¬ 
flares (with Swift) to highlight any possible change in the ac¬ 
cretion flow properties and/or geometry. In Figure we show 
the hardness ratio results for the Swift/XRT observations: in 
both outbursts the HR shows variability which seems corre¬ 
lated with the variations in luminosity, i.e. the source becomes 
softer at lower luminosities. 

Applying a standard Spearman’s rank test we find an cor¬ 
relation coefficient of 0.93 and 0.72 (for the 2008 and 2005 
outbursts respectively) but we note that the errorbars associ¬ 
ated with the colors and luminosities are significant and unac¬ 
counted for in this method. Using a composit e Monte Carlo 
analysis of Spearman’s rank test - as detailed in lCurranl (120141) 
and implemented in the code MCSpearman (Curran 2014) - 
we find a correlation coefficient of 0.71±0.21 and 0.66±0.18 
(for the 2008 and 2005, respectively) . This method uses a 
Monte Carlo re-sampling and perturbation analysis to return 
a robust correlation coefficient that accounts for both the un¬ 
certainty associated with the sample and the errors on indi¬ 
vidual data points. We therefore see a significant correlation 
between X-ray luminosity and HR, although we cannot draw 
any robust conclusion about the the lowest luminosities. 

In our fits to the SwiftlXKY spectra from the reflaring por¬ 
tion of the 2008 outburst (see Table 13 we used a spec¬ 
tral model consisting of some combination of a power law 



200 300 400 500 600 700 800 

Wavelength (nm) 



H I V b u w1 m2 w2 

Optical/UV Filter 


Figure 6. OpticalAJV flux vs. UVOT filter central frequency. In the top 
panel we show the flux Foptluvi'^) corresponding to one obseiwation taken 
during the main outburst. The dashed and solid lines represent the best fit 
to the data for F^p^iuy oc v“ (with a a free parameter) and F^pfUjy <=< v^, 
respectively. The results of the fit are reported in Table|4] 


(pegpwrlw), a blackbody (bbodyrad), and a disk blackbody 
(diskbb), modified by photoelectric absorption. We used the 
phab s absorption model with w ilm abundances ([Wilms et alJ 
|2QQQ|) and vern cross sections jVerner et^ll996b : the ab- 
sorption column was in all cases fixed at a value of A^h = 
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Figure 7. NIR/Optical/UV emission observed during the reflares (the sym¬ 
bols, with colors corresponding to different days, which are reported in Ta- 
ble[4) and and the predicted blackbody disk emission for the X-ray luminosity 
(continuous curves, with colors representing different days). Each curve cor¬ 
responds to a specific M inferred from the X-ray luminosity via the relation 
Lx GMMlR^ 



54755 54765 54775 53540 53570 53600 53630 


Time [MJD] Time [MJD] 


Figure 8. Hardness ratio (top panels; defined as the ratio between the count 
rate in the 2-10 keV and 0.3-2 keV band) of the 2008 (left panels) and 
2005 (right panels) reflares. On the bottom panels we report the ^wi/r/XRT 
lightcurves for comparison (we use only PC mode data and select only obser¬ 
vations that gave useful constraints on the hardness ratio). The HR and X-ray 
flux are correlated in both outburst: the spectrum becomes softer at lower 
luminosities. 

1.4 X 10^' Isee iPatruno et al.ll2009al) . For our spectral 
fits and count-rate conversions of the Swift data from the re¬ 
flaring portion of the 2005 outburst we used (in all cases) a 
model consisting of only an absorbed power law. The disk 
blackbody model is required only in the two brightest obser¬ 
vations of 2008. The choice of a disk blackbody is justified 
because such com ponent is also observed during the main 
outburst (see e.g.. IPatruno et al.l l2009al: iPanitto et al.l 120091: 
iKajava et al.ll201 ill . The results of our analysis are reported 
in Table |5] 

The power law index F shows little variations and all photon 
index values can be well fitted with a simple constant F = 
1.75 ±0.06 with a X^l^of — b.Sljl. However, given the 
poor photon counts and large error bars, the behavior at the 
lowest luminosities remains mostly unconstrained. 

We detect a blackbody with a small radius, which suggests 
we are observing the hot-spot or the neutrons star boundary 
layer rather than th e thermal emission from the inner accre- 
tion disk (see also IPatruno et al.ll200^ IPanitto et al.ll200^ 
for a similar analysis on XMM-Newton data taken during the 


main outburst). The two brightest reflare observations that 
need a multi-temperature disk blackbody (diskbb) still re¬ 
quire a rather cold inner disk temperature (0.1-0.2 keV). 

SAX J1808.4—3658 has remained in the hard state 
throughout its outburst history (from 199 8 until 2011) with 
very l ittle spectral variations (see e.g., Ivan Straaten et all 
I 2 OOI iH^an et al.1 l200i iBult & van der KlisI |20l4 . We 
caution that the minimum X-ray luminosity observable with 
RXTEfPCA is of the order of a few times lO^^ergs^^ 

4. DISCUSSION 

In this paper we have presented a multi-wavelength analysis 
on the accreting millisecond pulsar SAX J 1808.4—3658 that 
can help to address two specific (and related) issues: 

• what is the geometry of the accretion flow during re¬ 
flares? 


• what is the origin of reflares and how are they affected 
by the presence of a magnetosphere? 


4.1. Accretion Flow Geometry 

To understand how the accretion flow behaves during re¬ 
flares it is important to consider first whether the inner regions 
of the accretion flow are geometrically thin or thick and then 
how the flow interacts with the neutron star magnetosphere. 

Since SAX J1808.4—3658 has a measured dipolar mag- 
netic field (at the po les) of 2 x 10^ G (Idi Salvo et al.l 120081 : 
Hartman et al.ll2008l) . it is possible to estimate the location 
of its magnetospheric radius. We first define the magneto- 
spheric radius rm as the point where the magnetic field is 
strong enough to enforce co-rotation of gas in a thin Keple- 
rian disk: 


/ 77^2 


2 \ 1/5 


= 23 X lO^ri 
P* 


( 2 ) 


R. 


yio^Gj (lo^ 


1/5 


cm 

M 


6/5 


2xl0-3sy V 1.6 X 10-10 Mgyr-I 


-1/5 


cm. 


In this equation, /r = B^,Rf^ is the magnetic moment of the 
star, B* is the neutron star magnetic field at the poles, is 
the neutron star radius, Tj < 1 is a dimensionless parameter 
characterizing the strength of the disk/field coupling and re¬ 
flects the strength of the toroidal magnetic field induced by the 
relative rotation between the disk and dipolar magnetic field, 
(P*) is the star’s spin frequency (period) and M the mass 
accretion rate through the disc. If we use the known parame¬ 
ters of SAX J 1808.4—3658 and we assume that the observed 
X-ray luminosity is a good proxy for the mass accretion rate 
then at the lowest luminosities r„, ~ 160 km. This valuJ*^ 
is significantly larger than the co-rotation radius rco = 31 km 
and it is even larger than the light cylinder radius ric ~ 120 
km. The co-rotation radius is defined as that point in the disk 
where the Keplerian rotational velocity of the gas is equal to 
the rotational velocity of the magnetosphere and as soon as 


T he definition of used in this paper is th e same as inISpruit & Taaiil 
119931) whereas the classical expression of r^i IPringle & Reeslll972l) . ob- 
tained by equating the gas and magnetic pressures, gives a somewhat larger 
value of about 500 km, which would further strengthen the argument pre¬ 
sented here. 
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Table 5 

Spectral Fits to SAX J1808.4—3658 During the 2008 Reflaring State. 


DiskBB + BB + PL 

Observation 

DiskBB Temp. 

DiskBB Radius 

BB Temp. 

BB Radius 

r 

/^0.3—lOkeV.unabs. 

Z+dof 


[keV] 

[km] 

[keV] 

[km] 


[lO-’-^ erg/s] 


30034041 


Ho™ 

o.4857;;:iil? 

4 75+^'"^' 
^•'■^-0.61 

^•^^-0.29 

684+j^ 

41.2/50 

30034042 

U.UV/_oo23 

1304‘“ 

0.550«;«« 

^•ou_05g 

1 07 +O.II 

490f]j 

55.8/51 

BB + PL 

30034035 



- 

- 

2.00 (fixed) 

0.62+-- 

- 

30034036 



0.50(fixed) 

0 0O+0-27 

i.:5U_o 26 

10+11 

7.9/10 

30034037 



^•^^-0.06 

1 ^-T-I-O.SS 
^•^'-0.41 

1 00+0.23 
l-oo_0.24 

2o;i^ 

10.2/14 

30034038 



fl 4Q+0.08 
^•^^-0.06 

2 40+®-'^^ 
^•^^-0.53 

1 51+0-30 

122+1^ 

20.3/18 

30034039 



0.35(fixed) 

fl 62+^ "^^ 
^•^^-0.62 

^•°^-0.20 

11 7+^'^ 

3.5/5 

30034040 



- 


2.00 (fixed) 

i-+?5 

- 

30034043 



0 21C+0.03 
^•J4_oo2 

1 05+0-34 

1 70 + 0.14 
'"^-O.IS 

18.6™ 

36.4/39 

30034044 



0.25(fixed) 

1 ig+0.50 

1 75+^-29 
'^-0.43 

4-3;!1:’6 

- 


fm > rco a centrifugal barrier s ets in and matter is assumed 
to be ejected from the system (llllarionov & Sunvaevlll97-^ 
iRomanova et al.|[2Q04li . Therefore when > Tco one can¬ 
not use Eq. [3 since the M inferred from the luminosity might 
become a bad indicator of the mass flowing in the accretion 
flow (for example because of mass ejection). This means that 
Tm cannot be as large as 160 km, or at least, if it is, it can¬ 
not be inferred from Eq.|3 If the disk were truncated outside 
Tco it could also be truncated outside the light cylinder r\c. In 
this case the magnetosphere should be devoid of matter and 
the radio pulsar mechanism should turn on with a strong pul¬ 
sar wind preventing furth er accretion (see e .g., IStella et ^ 
11994 iBurderi et al.fcoOlh . This is what is currently thought 


sars recentlv discovered (lArchibald et al.l 

20091: iPanitto et al. 

2013 

; iPatruno et all 120141 IStacuers et al. 

120141 iBassa et al. 


2^ iRov et all 12014 and n the quiescence phase of 
SAX J1808.4-3658(lHomer et al.ll200U IBurderi et al.ll200l 
ICampana et all[2^ . The fact that in SAX J1808.4-3658 
the X-ray luminosity increases by three orders of magnitude 
right after reaching the luminosity minima on a very fast 
timescale of 1-2 days (see Eigure |2] and O suggests that 
the radio pulsar mechanism does not turn on, although a very 
rapid switch cannot be excluded at the moment. Eurthermore, 
pulsations at high luminosity (in SA X J1808.4—3658) and 
low luminosity (in PSR J1023-I-003 8. lArchibald et ^I20l4 
and XSS J12270-4859. iPapitto et alJl20l4 argue that accre¬ 
tion onto a neutron star can take place across a large range of 
accretion rates and therefore accretion onto the neutron star 
does not need to necessarily stop during the lowest luminos¬ 
ity phases of the reflares. 

In the remainder of this paper we will assume that the radio¬ 
pulsar mechanism does not turn on during the reflares (with 
the aforementioned caveats) and we will consider three possi¬ 
ble scenarios, sketched in Eigure |9] to understand what hap¬ 
pens to the accretion flow during reflares: 


A. : The disk is truncated well beyond the co-rotation radius 

and a strong propeller is ejecting the majority of matter 
flowing in the inner regions of the disk. 

B. : The inner disk is close to the co-rotation radius through¬ 

out the reflaring phase. 




Figure 9. Different possible accretion states in SAX J 1808.4—3658 . The 
vertical dashed line refers to the co-rotation radius rco- Panel A.: strong pro¬ 
peller configuration, when the disk is truncated well beyond rco and a very 
large fraction of gas is expelled from the system in an outflow due to the mag¬ 
netic centrifugal barrier. Episodic accretion is still expected to take place due 
to magnetic diffusivity, according to MHD simulations (see main text). Panel 
B. \ Trapped accretion disk scenario. The trapped disk forms when the mag- 
netospheric radius reaches rco and the disk remains truncated at that radius 
despite variations in the mass accretion rate. The accretion flow is modi¬ 
fied in the inner disk regions due to the centiifugal barrier, forming a “dead 
disk” configuration (purple segment of the disk). The outer disk regions ai‘e 
unaffected by the magnetic barrier and the disk is still a standard Shakura 
& Sunyaev disk (dark green segment). Panel C: Truncated/radiatively in¬ 
efficient accretion flow configuration. The inner disk region is set by the a 
mechanism other than the magnetosphere (e.g., evaporation). The resulting 
(geometrically thick) hot flow will still interact with the magnetosphere pos¬ 
sibly generating accretion and/or outflows. 
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Figure 10. Surface density profile times kinetic viscosity for an accretion 
disk for two different mass accretion rates. The red lines are the total ac¬ 
cretion disk density profiles, the dashed lines are two standard Shakura & 
Sunyaev accretion disk profiles and the dot-dashed lines refer to a dead disk 
density profile. At lO“^^M 0 yr“^ the accretion disk is dominated by the 
dead disk solution until r/rco ~ 10 and by the standard Shakura & Sunyaev 
profile beyond that point. At larger accretion rates, however, only the re¬ 
gions closest to the co-rotation radius are of the dead disk type and most of 
the accretion disk density profile is well described by a standard Shakura & 
Sunyaev solution. 


disk radius is set by the magnetospher^H If the inner¬ 
most region of the flow is not a thin disc, but instead 
forms a radiatively-inefficient accretion flo w (such as an 
advection dominated accretion flow, ADAF; |NaraYan _&JYi| 


1995:lNaravan & McClintock 20081 lAbramowicz et alJll995L 


20001) ■ then the situation might be somewhat different (sce¬ 
nario C). In most radiatively inefficient accretion flow (RIAF) 
models, the accretion flow becomes optically thin b ut geo¬ 
metri cally thick and somewhat sub-Keplerian (see e.g.. lSt)ruTd 
I 2 OOOI) . No detailed study of the interaction between a 
thick disk and a magnetosphere has been undertaken, and 
there is considerably uncertainty about the flow geome¬ 
try, location o f the magnetospheric radius a nd propeller ef- 
ficiency ( e.g., iMenou & McClintockI 1200 IL iDairOsso et alJ 
120141 and ID’ Angelo et al .1120151) . lit this case we expect a 
scenario somewhere between a thin-disk/magnetosphere in- 
terac tion and quasi-sphe rical accretion from a stellar wind 
(e.g. JShakura et al.l2013h . although there are many uncertain¬ 
ties about the properties of such a flow (e.g., the thickness of 
the hot flow, the inner ac cretion disk location, the pres ence 
of a jet and outflows; see lNarayan & McClintockil2008l for a 
review). To assess the validity of this scenario, further inves¬ 
tigations and theoretical modeling are required. 


C.: The thin disk is truncated far away from the neutron star 
magnetosphere which is interacting with a geometri¬ 
cally thick and radiatively inefficient hot flow. 

The first scenario (A) requires that a centrifugal bar¬ 
rier develops as soon as > r^o- To enter a strong 
propeller regime with matter expelled from the system, 
Tm needs to be larger than rco by at least a factor 
«1.3 so that the gas can g ain enough kinetic e nergy 
to reach the escap e velocity d Spruit & TaamI 11 9931). If 
a strong propeller (lUstyug ova et akl 120061 iRomanova et al.l 
120051 llllarionov & Sunvaevlll975T) is currently operating in 
the system then the very-low luminosities we observe are not 
due to an extremely low mass flow in the disk, but instead re¬ 
flect the fact that only a tiny (poorly-constrained) fraction of 
the mass actually falls onto the neutron star surface, generat¬ 
ing X-rays. In this case, the flow itself must likely be radia¬ 
tively inefficient since otherwise the much-higher accretion 
rate in the disc will likely dominate the X-ray emission. In 
this specific case the term radiatively inefficient is not used to 
indicate the presence of a geometrically thick radiatively inef¬ 
ficient flow (like an advection dominated accretion flow) but 
it is used to describe the fact that the majority of the accretion 
energy is not emitted as radiation. 

If Tm remains too close to rco so that matter cannot be 
ejected from the system, then an accretion disk different than 
the standard geometrically thin and optically thick Shakura- 
Sunyaev type needs to develop (scenario B). Since the cen¬ 
trifugal barrier is not sufficient to expel the gas, the angular 
momentum is transferred from the magnetosphere to the mat¬ 
ter that now remains in the inner disk regions. The gas piles 
up and changes the radial density and temperature profile of 
that region (see e.g.. Figure fTOll. This segment of the disk is 
then dominated by a trapped/dead disk solution rather than 
a stan dard Shakura and Sunyaev one (ISiuniaev & Shakti^ 
Il977h . In this case the X-ray luminosity is low because only a 
small fraction of matter might leak from the inner (dead) disk 
towar ds the neutron star surface (ID’ Angelo & Spruidl20ll 

IMl) . 

Both these scenarios require that the inner accretion 


4.1.1. Observational Constraints on Accretion Flow Geometry 

Having set the three different configurations for the geom¬ 
etry of the accretion flow and disk/magnetosphere interaction 
during the reflares, we now discuss how the new observations 
presented in this work and the past observational evidence 
collected from the literature can help to constrain the differ¬ 
ent scenarios. The three possibilities discussed above have a 
number of specific testable features that can be translated into 
the following questions: 

• is there a geometrically thin disk extending all the way 
down (or close to) the co-rotation radius? Or is the in¬ 
ner disk radius truncated much further away from the 
neutron star? 

• does the disk/magnetosphere interaction generate a 
strong propeller with outflows? 

Constraints from X-Ray Spectral Analysis. 

Our spectral analysis of the reflares indicates that 
SAX J1808.4—3658 has an overall stable spectrum de¬ 
spite the large variations in luminosity. The Swift/XRT 
spectra remain hard at all luminosities and all spectra are con¬ 
sistent with a constant value of F ~ 1.7. We do find correlated 
variations between the Swift hard color and X-ray luminosity 
that suggests that some (small) changes in the system are 
occurring in response to the X-ray luminosity variations. 
The RXTE colors further show that SAX J1808.4—3658 
is observed in the hard (island) state throughout its eight 
outbu r sts (see e.g., lin ’t Z and et al.l 119981: Iv an Straaten et al.l 
120051: iHartman et akl l2008t iBult & van der KlisI 120141 and 
Section I3.41 i. The rather stable X-ray spectrum suggests 
that the production of X-rays is regulated by some stable 
mechanism despite large variations in the accretion flow. 

During the brightest portions of the reflares we detect an 
accretion disk in addition to a blackbody and a power-law. A 

Tlie region over wtiicti ttie dislc and magnetic fieid interact will of course 
have a certain radial extent. However, both simulations and t heoretical argu¬ 
ments predict this transition is sud den and (see e.g.,ILovelace et al.l 

119951 and an extended discussion in ID’Angelo & Spruitll2Ql(]n . 
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very similar spectrum was seen also during the main 2008 out¬ 
burst with XMM-Newton ([Patruno et alj| 2009at iPanitto et alJ 
[^(MlCackett et al.ll2009HKaiava et al.ll207llK 

An important caveat is that even if the spectra of the reflares 
are similar to the main outburst spectra, this does not neces¬ 
sarily mean that the accretion flow has the same config ura¬ 
tion. Indeed the spec tr a analyzed by Kaiaya et al.l (1201 lb and 
iPatruno et al.l (l2009al) . [Panitto et al.l ( 2009^ and ICackett et aT] 
(120091) refer to luminosities about one order of magnitude 
brighter than the reflares. 

Constraints from X-Ray Aperiodic Variability. 

A truncation radius close to the co-rotation one is also 
the preferred scenario to explain the strong 1 Hz modula¬ 
tion observed during the 2000 , 2002 and 2005 reflares of 
SAX J1808.4-36 58 (IPatruno et al.1 l2009bt IWiinandsI [2001 
Ivan der Klisll2000l) . Recently, a similar modulation (1-5 Hz) 
has been observed also during the main outbu rst of the 2008 
and 2011 outbursts (iBult & van der Kl'isll2014l). In both cases 
an instability arising from a trappe d disk dP’Angelo & Spruill 
1201 Ol 1201 21 : ISpruit & TaamI 119931) has been proposed to be 
at the origin of this phenomenon. Such a disk model has 
also been discussed to explain a strikingly similar 1 Hz 
modulation observed i n another accreting millisec ond pulsar 
(NGC 6440 X-2; see IPatruno & D’ AngelollMl . It is not 
possible to exclude that this variability mig ht be related to 
propeller driven episodic accretion (e.g., iLii et alJ IMl) 
which, however, needs further investigation. 


Constraints from Coherent Pulsations. 

Even if there are indications for the presence of a thin disk 
extending all the way down to the co-rotation radius, we still 
do not know whether we are in the presence of a propeller (A) 
or a trapped disk (B). Important indications can be given by 
the observations of coherent pulsations. 

The detection of pulsations during reflares (up to the 
sensitivity limit of RXTE of about a few lO^^ergs^*, see 
IPatruno et al.ll 200 ^ and the presence of a blackbody in the 
spectrum (down to luminosities of lO^^ergs^*, see Table |5]l 
does indicate that some gas keeps falling onto the neutron star 
surface during most of the reflares. Therefore accretion on the 
neutron star surface is not (completely) inhibited during re¬ 
flares, most likely at least not down to ^10^^ ergs^^ This is 
compatible with the recent findings of accretion powered pul¬ 
sations in two quiescent LMXB, where pulsations are seen at 
Lx — lO^^ergs^* (lArchibald et alJl2014l : lPapitto et al.ll2014^ ■ 
so that accretion at these luminosities is not an unprecedented 
occurrence. The presence of pulsations and a blackbody at 
low luminosities imply that - if a strong propeller is currently 
operating (rather than a “trapped disk”) - sustained accretion 
onto the neutron star surface is still ongoing. Numerical MHD 
simulations show that during a strong propeller only a small 
portion of the mass flow (of the order of few percent, although 
the exact value strongly depends on the poorly constrained 
coupling between the magnetic held and the plasma) actu¬ 
ally reaches the neu t ron star surface (e.g., v i a episodic accre¬ 
tion, Lii et al.l 120141; lUstvugova et alj|2006t I Romanova et alJ 
l2005l) creating a boundary layer/hot-spot that gives rise to 
most of the observed X-ray luminosity. Therefore if a strong 
propeller is currently operating in SAX J1808.4—3658 we 
should expect a very large mass outflow. 

Such large outflow of mass might have some further ob¬ 
servational consequences. Indeed if the amount of mate¬ 
rial reaching the surface is truly only a very small fraction 


of the total amount of matter flung off the disk, then we 
should still expect a rather strong spin down of the neutron 
star. So far, SAX J1808.4—3658 has not been observed to 
spin up or down during its outbursts, with typical upper lim- 
its of Ivl < 2.5 X IQ-^^Hzs-^ (IHartman et alJl2(M l2(M 
IPatruno et al.ll2012h . a fact that again suggests that this object 
does indeed stay close to spin eq uilibrium (rm—rm) durin g 
most of its outbursts (see also e.g.. lHaskell & Patrunoll20ll]) . 
Although these upper limits refer to the entire outburst and not 
only to the reflaring portion, we note that at least in one case 
(the 2000 outburst) RXTE collected data exclusively during 
the reflares (due to solar constraints). In that case the upper 
limits on the spin frequency derivative were even more strin¬ 
gent: -1.1 X lO-'^^Hzs-' < V < 4.4 X IQ-'^Hzs-^ (95% 
confidence interval, see Table 4 iii THartman et alj|2008l) . 

Episodic accretion might possibly have a role to explain 
why we do not see a strong spin down despite the strong pro¬ 
peller since the angular momentum loss due to the centrifu¬ 
gal barrier might be (partially) compensated by the material 
torque accreting onto the neutron star surface. However, to 
determine the order of magnitude expected we can use the re¬ 
lation: 

4op = -nMej(GMnn)‘/2 

= -«(nn/rco)‘/^Mej(GMrco)‘/2, (3) 


where rin is the inner disk radius and the dimensionless 
torque n is zero for rjn = rco and of order unity for rjn > 
1.1 rco (lEksi et al.l l2005l) . Eor rin « 1-2 rco and Mej ~ 
lO^*’'M 0 yr^^ we obtain a spin down of the order of a few 
times 10^'^Hzs^*. Such values are of the same order of mag- 
nitu de as the upper limits placed in mo st of the outbursts (se e 
e.g.. IHartman et aO2009l and Table 4 in iHartman et alJl2008l) . 
The choice of Mej ~ 10^^** M© yr^* derives from the fact that 
this is the approximate value required to produce X-ray lumi¬ 
nosities of a few times 10 ^^ ergs^*, which is the typical value 
observed at the peak of the reflares. Since the material falling 
on the neutron star has to be at least 10 ^"^ M 0 yr^*, then the 
expected spin-down above is the minimum that we can expect 
in this scenario. The constraints on the observations (at least 
for the 2000 reflares) suggest that any spin down present in 
SAX J1808.4—3658 must be smaller than this minimal ex¬ 
pected V value. However, given the uncertainties in the model 
and the close value between the expected v and the observed 
upper limits, this result should be taken more as a possible 
indication against the strong propeller rather than conclusive 
evidence. 

The lack of a detected spin-down can also be important to 
place constraints on the trapped disk scenario (B). The sign 
of the torque depends on whether the inner edge of the disk is 
inside or outside the co-rotation region, and the amplitude is 
roughly determined by the “critical accretion rate” - Me, the 
accretion rate at which th e inner edge of the disk is equal to 
the co-rotation radius. In ID’Angelo & Spruitll2012l this was 
estimated as: 


Mc = 


4n*r5’ 


(4) 


and corresponds to Me 5 x 10 yr ^ for 

SAXJ1808.4-3658 (assuming a magnetic held of lO^G). 
The amplitude of the spin-up/spin-down torque is then 
roughly: 

./c=Mc(GM*rc)‘/^ 


(5) 
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which (assuming a moment of inertia of cm^) corre¬ 

sponds to a spin change of Av 10^'^Hz. For a mean out¬ 
burst accretion rate of M ^ 2 x Mq yr^' this corre¬ 

sponds to a spin up rate of ~ (0.3-1) x 10^'^Hzs^^ For 
the long-term average accretion rate of M ^ 5 x 10^*^ Mq 
yr^* (i.e., averaged over an entire outburst/quiescence cy¬ 
cle), the spin-up from accretion is nearly balanced by spin- 
down from the disk-magnetosphere coupling outside rc, and 
more detailed calculations predict a net spin down rate of (1- 
3) X 10~^^F[zs~* (see also Figure 2 of [D’Angelo & Spruill 
l2012h . Such values are not much different than those obtained 
for the propeller scenario (A) and indeed similar conclusions 
can be drawn: the expected spin up/down from a trapped disk 
is too close to the existing measured upper limits to place ro¬ 
bust constraints on the operating mechanism. 

Finally, it is worth noticing that the lack of a spin-down 
during the outbursts might still be compatible with the 
scenario C, where a thick hot flow is present up to a very 
large distance from the neutron star. In this case it might be 
the sub-Keplerian nature of the hot flow itself that results in a 
very small amount of angular momentum transferred towards 
the neutron star and thus in a lack of detected accretion 
torques. In Table 6 we summarize which of the observed 
properties can be explained by each of the models considered 
here. 


4.2. The Origin of Reflares 

Reflares are an important phenomenon to study because 
they are assumed to be generated by the ionization instability 
that drives dwarf novae and X-ray binary outbursts. There¬ 
fore understanding why reflares are sometimes observed can 
give important information about how the ionization instabil¬ 
ity operates. Even more importantly, their observation can 
help testing the assumption that they are indeed generated by 
the ionization instability. 

Reflares are currently considered a problem for the disk 
instability model since they require that a large amount of 
matter is retained in t h e disk a the end of the main out¬ 
burst (lArai et alJ 120091: iDubus et'aD 1200 Ih . The currently 
best explanation (which qualitatively agrees with the observed 
timescales and luminosity fluctuations) is that reflares are es¬ 
sentially ‘mini-outbursts’. A small change in disk density at 
the end of an outburst increases the outer disk temperature 
enough to partially ionize hydrogen, which then leads to a 
rapid rise in the accretion rate into the inner disk and a re¬ 
brightening. The exact trigger of the reflare is uncertain - 
they a ppear spontaneously in the simulations of [Dubus et alJ 
(120011) . although they do not resembl e the o bserved reflares 
and they are seen by iHameurv et alJ (l2000l) where reflares 
are caused by an increased irradiation of the donor star that 
causes a superoutburst (so called because their duration is 
much larger than that of normal outbursts) followed by re¬ 
flares. The donor star in SAX J 1808.4—3658 is observed to 
be strongly irradiat e d during quiescence (Homer et al.l120011 : 
iBurderi et alJ 120031 : ICampana et alJ 12004 ) and indeed there 
are suggestions that it is losing a large amount of m ass 
(Idi Salvo et al.l 1200^ see however iPatruno et al.l 120121 and 
iHartman et al.l [200^ for criticisms of this strong mass loss 
scenario in SAX J1808.4—3658). Another suggestion for 
the re flare trigger comes from the mass reservoir model of 
lOsaki et aP (120011) . where reflares are triggered also after su¬ 
peroutbursts as long as the effective viscosity of the disk 


(parametrized by the a parameter) remains large through the 
entire sequence of reflares. 

So far little work has been done on reflares occurring in X- 
ray binaries and most theoretical work has focused on explain 
the reflaring phenomenon observed in the dwarf novae pop¬ 
ulation. However, regardless of the trigger (which might be 
of a different nature than outlined above), the reflares should 
proceed in the same way: at the end of an outburst the temper¬ 
ature at some location in the outer disk falls below the hydro¬ 
gen ionization temperature, causing hydrogen to recombine 
and the accretion rate to fall dramatically as a ‘cooling wave’ 
travels inward. For (possibly) one of the reasons outlined by 
the irradiation or mass reservoir model, the density in the disk 
subsequently rises again enough to reach the hydrogen ioniza¬ 
tion temperature and again send a new ‘heating’ wave through 
the disk, which causes the accretion rate to rise. Since there 
is not a large amount of gas left in the accretion disk after the 
main outburst, the accretion rate again decreases, resetting the 
entire process. 

In WZ-Sae type systems, the reflares are observed after the 
occurrence of a superoutburst, so called because its duration 
is much larger than that of normal outbursts. In these systems, 
it was proposed that a resonance occurs when the disk radius 
exceeds the threshold of 0.46a where a is the orbital semi¬ 
major ax is. When th i s happ ens superhumps are expected to 
develop (lOsakill 198911199^ . Black hole candidates in X-ray 
binaries displayi ng superhurnps do not show outbursts of var¬ 
ied amplitudes (iMaccaronel 120141) and SAX J1808.4—3658 
does follow the same behaviour. It is possible that both 
SAX J1808.4—3658 and those black hole candidates have 
shown only superoutbursts. iMaccarone & Patrunol (1201 3h 
suggested indeed that we do observe only super-outbursts 
from some LMXBs and superhumps should be detected if the 
reflares in SAX J1808.4—3658 are a variant of those seen in 
d warf novae. _ 

lElebert et akl (l2009h reported weak evidence for the pres¬ 
ence of a superhump at optical wavelengths in SAX J1808.4- 
3658, which, if confirmed, might be an exciting observational 
diagnostic and might indicate the presence of an eccentric 
disk. However, the candidate superhump has been observed 
so far only once during the main 2008 outburst and it attends 
further confirmation. 

Regardless of the presence of superoutbursts, this basic 
hot/cold wave scenario can plausibly account for the reflares 
in SAX J 1808.4-3658. If the instability originates in the outer 
part of the disk, the reflare rise time can be associated with the 
time it takes for the heating front to propagate from the outer 
edge of the disk to the central regions, which will happen on 
a thermal timescale: 

fr -- rd/acs (6) 

where rd is the location of the outer edge of the disk, a = 0.1 
is the Shakura-Sunyaev viscosity parameter, and Cs is the 
sound speed at the inner edge of the heating front, where 
T ~ 6500 K. Eor a 2-day rise time, rd ^ lO'** cm, which 
is comparable to the circularization radius of the disk r^ ~ 
0.46a ^ 3 X 10*** cm (assumed to be the outer disk edge; a 
is the projected semi-major axis of the binary, for the orbital 
period 7249s and assumed NS mass of 1 .4Mq and donor com¬ 
panion of 0.01 Mq). 

4.2.1. Constraints from Multiwavelength Photometry 
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Table 6 

Models for the Geometric Configuration of the Accretion Flow 


Model 

X-Ray Spectra 

Spin-Down 

1 Hz QPO 

Weak Irradiation 

X-Ray Delay 

NIR/Optical/UV Excess 

Strong Propeller (A.) 

P 

P 

P 

? 

Y 

Y 

Trapped/Dead Disk (B.) 

Y 

P 

P 

? 

Y 

Y 

Thick Hot Flow (C.) 

N 

P 

N 

? 

Y 

Y 


The table compares observed properties of the to the various models discussed in §4. Y/N (yes/no) indicates that the model can/cannot 
explain the property in SAX J 1808.4—3658. The symbol P (possible) indicates that the model might be able to explain the property if certain 
conditions are met. The symbol ? indicates that the model makes no specific prediction for that property, and that further studies are required. 


During the reflaring phase of the 2008 (and to some extent 
also 2005) outburst, we can place three constraints on the na¬ 
ture of the accretion disk in SAX J1808.4—3658. The first 
constraint comes form the observation of wea k irr adiation as 
witnessed by the small value of j3 (see Section lT^ . This usu¬ 
ally suggests that the accretion disk is heated mostly by vis¬ 
cous dissipation rather than irradiation (see e.g.. lRussell et al.l 
1 20061) . The reason why irradiation is not the main contributor 
to the optical/UV flux emission remains an open question. 

The second constraint comes from the observation of some 
anti-correlations between X-ray and optical/NIR emission 
(e.g., Figure |2] on MJD 54746-54747 and MJD 54758 and 
several data point in Figure |3ll. As discussed in Section [TTI 
there is a tentative X-ray delay of about 1.5 days with re¬ 
spect to the NIR/optical emissioiQ The X-ray delay is a phe¬ 
nomenon that has been already observed in other LMXBs al¬ 
though so far only in systems containing black holes. In some 
systems (two of which are persistent sources) the opticaFX- 
ray cross-correlations sho ws evidence for X-ray delays on 
the order of ~5-20 days. iBrocksopp et al.l (120011) found an 
X-ray lag of 5-10 days from long-term optical and X-ray 
light-curves of the b lack hole binary L MC X-3 (a persis¬ 
tent source; see als o ISteiner et al. J20l4 who measured this 
delay as 2 weeks). iHoman et al.l ( 200*^ found an X-ray lag 
of 15-20 days compared to near-IR variations in the black 
hole LMXB GX 339-4 (a transient source). An X-ray lag of 
2-14 d ays was observed in 4U 1957-1-11 (persistent) behind 
optical (iRussell et al.ll2010l) . whereas a 3.5 d X-ray time lag 
compared to UV was found in Sw ift J1910.2-0546 (a tran¬ 
sient source; iNakahira et al.l 1201^ . Such delays have been 
interpreted with a variety of models that propose that the 
NIR/opticaFUV emission originates from a jet, from the outer 
disk or even from the irradiated stellar surface. 

In the case of SAX J1808.4—3658 we do not have suffi¬ 
cient observational evidence to constrain the presence of a jet 
during the reflares. Since the jet model (or an ADAF with 
an outflow) cannot be constrained we will not consider it any 
further, although we stress that there is still no proof that can 
firmly exclude such scenario. 

If we consider a disk origin of the emission then we can 
assume that the NIR/optical/UV emission is coming from the 
outermost disk regions. In this case by looking at Figure|2]it 
is clear that this emission is over-luminous with respect to the 
X-rays luminosity. It is possible that, when the heating front 
moves inward and reaches the inner disk regions, a larger flow 
of X-rays is produced which would correspond to the peak of 
the reflares. However, the peaks are still too dim for the opti¬ 
cal luminosity observed and therefore this scenario is unable 

We also note tha t during some o bservations of the 2000 outburst 
(see lWachter et alll2000l : rWiinandsir200fD the optical luminosity also did not 
con'elate with the X-ray luminosity. 


to explain the observed behaviour. 

On the other hand, if an outflow is present and a signifi¬ 
cant portion of mass is lost from the system then a strongly 
over-luminous NIR/optical/UV emission is expected. Indeed 
if the material generating radiation is lost in the disk before 
reaching the innermost regions then the X-rays will not be 
produced. This can explain the possible X-ray delay, as well 
as our third important constraint: any multicolor-disk model 
of a viscously heated accretion disk seems to modestly under¬ 
predict the optical luminosi ty ob served for the measured X- 
ray luminosity (see Section (33]). Since the optical emission 
is only over-predicted by a factor of a few with respect to 
the X-ray flux (assuming a constant accretion rate through 
the disk) our results also exclude a very strong outflow in 
the disk. Such an outflow would expel the majority of gas 
before it radiated energy in X-rays, which would make the 
discrepancy between optical and X-rays much larger. How¬ 
ever, it is still possible that some (not too strong) outflow 
of material is present in the system, which could be driven 
either by a magnetosp heric RIAF (like the strong propeller; 
lUstvugQv a et alJl 2006t) orbv an accr etion driven RIAF (like a 
windi lBlandford & BegelmMlll999h . is present in the system. 

Finally, almost identical conclusions can be reached for the 
scenario B (trapped disk). In Figure [TO] we compare the sur¬ 
face density profiles (times the kinetic viscosity) of a stan¬ 
dard accretion disk and a trapped/dead disk. The extent of 
the trapped disk region could be modest which means that 
the largest majority of the accretion disk would still look 
and behave as a standard Shakura & Sunyaev accretion disk. 
This means that at large accretion rates the total disk profile 
would be observationally indistinguishable (in luminosity and 
spectral properties) from a Shakura & Sunyaev disk, despite 
the presence of a small segment where the dead disk is still 
present. Therefore identical considerations discussed above 
apply to this scenario too. We caution that the extent of the 
region where the dead disk dominates strongly depends on 
the efficiency of the angular momentum extraction at the in¬ 
ner edge of the disk, a quantity that is currently poorly known. 

Interpreting the reflares in terms of a specific physi¬ 
cal model is hampered by the uncertain origin of the 
UV/opticaFNIR emission. Assuming that it is dominated by a 
single physical component, the simplest explanation would be 
that it originates in the outer parts of the accretion disk, either 
directly from viscous dissipation or as reprocessed emission 
from irradiatio n from th e inne r X-ray component. However, 
as discussed in 13.31 and 14.1.11 we consider it unlikely to be 
reprocessed emission and its origin remains still to be identi¬ 
fied. 

5. CONCLUSIONS 

In this paper we have analyzed the 2008 and 2005 
reflaring phase of the accreting millisecond pulsar 
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SAX J1808.4—3658. We find that the accretion flow 
geometry can be substantially constrained by the occurrence 
of reflares. We find that a disk extending down to the 
co-rotation radius and truncated by the neutron star magneto¬ 
sphere provides the most likely explanation for the accretion 
flow. We have found that the NIR/optical/UV emission 
during the reflares is over-luminous and cannot be explained 
by an irradiated accretion disk model if the mass accretion 
rate used is inferred from the observed X-ray luminosity. We 
suggest that most of the NIR/optical/UV emission comes 
from the outer disk and it is produced by viscous dissipation. 
The mass accretion rate in the outer disk regions needs to 
be substantially higher than the one that reaches the inner 
disk regions. We propose that either a propeller with a 
large amount of matter being expelled from the system or 
a trapped (dead) disk truncated at the co-rotation radius are 
both plausible explanations for all of these features. A disk 
truncated (by e.g., evaporation) far away from the neutron 
stars (at hundreds or even thousands gravitational radii) is 
instead not supported by the current observational evidence. 
Finally, we find that a hot/cold wave propagation model for 
the reflarings is compatible with the observed timescales 
although the mechanism that triggers such phenomenon is 
still to be identified. 
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